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Abstract 

A variety of resting state neuroimaging data tend to exhibit fractal 
behavior where their power spectrums follow power-law scaling. Resting 
state functional connectivity is significantly influenced by fractal behav- 
ior which may not directly originate from neuronal population activities 
of the brain. To describe the fractal behavior, we adopted the fraction- 
ally integrated process (FIP) model instead of the fractional Gaussian 
noise (FGN) since the FIP model covers more general aspects of fractal- 
ity than the FGN model. This model provides a theoretical basis for the 
dependence of resting state functional connectivity on fractal behavior. 
Inspired by this idea, we introduce a novel concept called the nonfractal 
connectivity which is defined as the correlation of short memory indepen- 
dent of fractal behavior, and compared it with the fractal connectivity 
which is an asymptotic wavelet correlation. We propose several wavelet- 
based estimators of fractal connectivity and nonfractal connectivity for a 
multivariate fractionally integrated noise (mFIN). These estimators were 
evaluated through simulation studies and applied to the analyses of resting 
state fMRI data of the rat brain. 

1 Introduction 

The dynamics of endogenous neuronal activities has been an important issue 
in neuroscience since it is supposed to take control of most neuronal activities 
arising in the brain 1 1 . The huge default-mode functional network of the brain 
has been usually investigated through resting state neuroimaging data such 
as electroencephalography (EEG) and functional magnetic resonance imaging 
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(£MRI) [2jj5]. One of the major goals in resting state neuroimaging research is 
the reliable inference of functional dynamics of spontaneous neuronal population 
activities from resting state neuroimaging data. However, it is not straightfor- 
ward since resting state signals may be significantly affected by non-neuronal 
physiological factors. On the other hand, the response to stimulation in task- 
based experimental paradigm is prominently correlated with brain activities 
either directly or indirectly. 

One of the non-neuronal obstacles in resting state neuroimaging studies is 
the fractal behavior (or long-range dependence) where the power spectrum tends 
to exhibit 1/ f a power law scaling across low frequencies. This phenomenon has 
been observed through a number of resting state neuroimaging studies [6jj9]. 
As such a phenomenon has been ubiquitously observed in nature, the fractal 
behavior in neuroimaging data may also arise from various mediators such as 
respiration [lOj[ll], cardiac fluctuations [l2], system noise, hemodynamics (in 
the case of fMRI) as well as neuronal activities [l3}{l5] . 

The classical model of fractal behavior or long memory in baseline neu- 
roimaging signals has been the fractional Gaussian noise (FGN) which is defined 
as an increment process of fractional Brownian motion (FBM) and completely 
characterized by both Hurst exponent and variance |16|. The FGN model has 
been adopted to various fractal-based analyses of fMRI data for a decade so as 
to account for scale- free dynamics of neuroimaging signals [l7]j2l] . 

However, there is a controversy about whether the FGN is the most ap- 
propriate model for resting state neuroimaging signals among a variety of long 
memory models. While the FGN model is defined just with two parameters 
under mathematically strict conditions of self-similarity, a neuroimaging signal 
is produced from a nonlinear biological system which is controlled by numer- 
ous hidden parameters. In this reason, the fractionally integrated process (FIP) 
model, based on the concept of fractional differencing [22], is worth considera- 
tion as an alternative to the FGN model since it embraces diverse types of long 
memory. Indeed, the FGN is regarded as a special type of the FIP model which 
is more extensive than FGN. 

In this paper, we adopted the fractionally integrated process (FIP) model 
to effectively describe the fractal behavior of neuroimaging signals. In the FIP 
model, a neuroimaging signal is represented as the output of a long memory 
(LM) filter whose input is a nonfractal signal (sometimes called short memory 
as a notion corresponding to long memory). In other words, a nonfractal signal 
is transformed into a neuroimaging signal with fractal behavior through long 
memory filtering as shown in Fig. [T] which indicates that the fractal behavior 
is attributed not to the nonfractal input but to the LM filter. The influence 
of several factors on the fractal behavior can be well aggregated in terms of a 
sequence of long memory filters. 

The FIP model sheds light on the influence of fractal behavior on functional 
connectivity. The correlation of resting state neuroimaging signals between 
two brain regions may significantly differ from that of the nonfractal input 
signals according to the difference of memory parameters. Hence, the ordinary 
correlation of resting state neuroimaging data may not well reflect the functional 
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Figure 1: The FIP model representation of resting state neuroimaing signals. 



dynamics of spontaneous neuronal activities due to the fractal behavior. 

The theoretical expectation that functional connectivity may be influenced 
by fractal behavior leads us to take into account the correlation of fractal-free 
input signals as a novel concept of resting state functional connectivity while the 
Pearson correlation of neuroimaging signals has been the most popular defini- 
tion of functional connectivity. This particular correlation, which is independent 
of fractal behavior, is called the nonfractal connectivity. Its mathematical de- 
scription is provided in section [3] The nonfractal connectivity is not exactly 
identical to the correlation of spontaneous neuronal population activities due to 
the nonlinearity of neurophysiological systems. However, it may give us better 
information on correlation structure of spontaneous neuronal populations than 
ordinary correlation of neuroimaging data since it eliminates the distortion of 
functional connectivity due to fractal behavior. 

The nonfractal connectivity is comparable to fractal connectivity which was 
first proposed in [23 as the asymptotic value of wavelet correlations over low 
frequency scales. The wavelet correlations of two long memory processes are 
converged on a specific value which is determined by the correlation of short 
memory parts as well as memory parameters. We show the theoretical relevance 
of nonfractal connectivity to fractal connectivity. 

In this paper, we propose three wavelet-based approaches to estimating both 
nonfractal connectivity and fractal connectivity: (1) the SDF (Spectral density 
function)-based method, (2) the covariance-based method, and (3) the linearity- 
based method. As prerequisite to estimating these connectivities, memory pa- 
rameters should be estimated a priori. We tested two wavelet-based univariate 
estimators of memory parameter: the wavelet-based least-mean-squares (LMS) 
method and the wavelet-based maximum likelihood (ML) method. The perfor- 
mance of all proposed estimators was verified through simulation studies. We 
also show an example of applying these estimators to resting state fMRI data 
taken from anesthetized rat brains to estimate nonfractal connectivity. 

This paper is organized as follows. In section |2j the FIP model of resting 
state neuroimaging signals is briefly introduced, and the concepts of nonfractal 
connectivity and fractal connectivity are described in section [3j The wavelet- 
based estimators of both connectivities are proposed in section [4j The results 
of simulation studies and experiments in fMRI data are provided in section [5] 
and [6l 
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2 Long memory model 



Fractal properties of a time series can be modeled as long memory models such as 
FGN, ARFIMA, and GARCH processes. In this section, we especially introduce 
both univariate and multivariate fractionally integrated processes (FIP) that 
encompass several classes of long memory such as fractionally integrated noise 
(FIN), FGN, and ARFIMA 



22 



2.1 Univariate case 

Let x(t) be a real- valued discrete process of length N given by 

u(t) = (1-L) d x(t) 



(1) 



where d G M, L denotes the back-shift operator, and u(t) (called short memory) 
is a stationary process whose spectral density S u (f) is a non-negative symmetric 
function bounded on (— 7r,7r) and bounded away from zero at the origin. x(t) 
can be represented as the convolution of u(t) with the long memory (LM) filter 
g(t) as follows 



x(t) = J29(r)u(t-r) 



r=0 



where 



g(t) 



dT(d + t) 



If 



r(d + i)r(t + i) 

-1/2 < d < 1/2, the spectral density of x(t) is given by 



S(f) 



1 



if\ 



-2d 



Su(f). 



(2) 



(3) 



(4) 



The fractal behavior is controlled by the memory parameter d. If0<<i<l/2, 
the process x(t) is said to be a stationary long memory process with memory 
parameter d while x(t) is nonstationary if <i > 0.5. If d = 0, the process becomes 
a white noise. 

2.2 Multivariate case 

The definition of the univariate long memory model can be extended to the 
multivariate case. Consider a real- valued q- vector process X(t) given by 



V o 







(1 - L) d « 



V x q (t) 



( Ui{t) 



\ US) 



(5) 



where U(t) = (C7i(t), U q (t)) is a multivariate stationary process whose spec- 
tral density S(/) = [Sm in (f)] is bounded on (— 7r,7r) and bounded away from 
zero at the origin. For —1/2 < dj~ < 1/2, the spectral density of U is given by 



S(/) = $(/)S„(/)$* (/) 



(6) 
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where 

/ (l-e^)" dl \ 

*(/)= • (7) 

\ (l-e if )- d « J 

In the case of < dk < 1/2 for 1 < k < q for 1 < k < q, X(t) is said to be 
a stationary long memory process with memory parameter d = (di, • • • , d q ). If 
U(t) is a vector ARM A process, X(£) becomes a multivariate ARFIMA process. 
On the other hand, if U(t) is a vector i.i.d. random variable, i.e. 

U(i) L ^ AT(0,£„), (8) 

X(t) becomes a multivariate fractionally integrated noise (mFIN). In this case, 
the cross-spectral density of x m (t) and x n (t) is given by 

S m ,n (/) = lm,n (l - e^) " dm (l - e"'/) (9) 
where jrn^n is identical to the (m, n)-th element of E n . 



3 Nonfractal and fractal connectivity 

As discussed in the section [T] the most popular definition of functional connec- 
tivity has been the Pearson correlation. The multivariate long memory model 
introduced in the section [2] additionally provides two novel definitions of resting 
state functional connectivity: fractal connectivity and nonfractal connectivity. 
While fractal connectivity is defined based on the asymptotics of wavelet cor- 
relation, nonfractal connectivity is defined based on the covariance of short 
memory. 



3.1 Nonfractal connectivity 

Let X(t) be an mFIN process with memory parameter d, and U(t) be a short 
memory function of X(t) given in (J5|. The nonfractal connectivity of x m (t) and 
x n (t) is defined as 

D m>n = 7m '" (10) 

V7m,m7n,n 

where 7 m , n denotes the covariance of u m (t) and u n (t); that is, 7 m?n := E [u m (l)u n (l)]. 



3.2 Fractal connectivity 

The variance of a discrete time series can be decomposed over several fre- 
quency bands (called scales) through the discrete wavelet transform (DWT). 
Let Wi(j, k) be the wavelet coefficient of the zth process Xi(t) at scale j and time 
point k. The wavelet covariance is defined as ^ m?n (j) := cor (W m (j, fc), W n (j, k)) 



5 



at scale j. Since the wavelet coefficients of an FIP at a scale j is covariance sta- 
tionary, Vm,n(j) is independent of time t. Let Hj(f) be the squared gain function 
of the wavelet filter such that 



f2^ for 27T/2J+ 1 < |/| < 2-k/V 
j[J ' ~ \ otherwise 



(11) 



Then, the wavelet covariance of x m (t) and x n (t) at scale j is related to the 
cross-spectral density |25 as follows 



VmAi) = 2tt ^ Hj(f)Sx(f)df. 

J —TV 

The wavelet correlation p m ,n(j) : = cor (Wm(j, /c), W n (j, k)) is given by 

^m,n(j) 



Pm,nU) 



(12) 



(13) 



Theorem 1 (Asymptotic wavelet covariance). Suppose thatX(t) is a multivari- 
ate FIN process which satisfies the condition (|8| . Then, the wavelet covariance 
of x m (t) and x n (t) at scale j is approximated by 



where 



^m,n{j) ~ ^rn,nfirn,n^^ n ^ as j 

]_ _ 2 d m+d n — l 



oo 



2 cos 



(§<4 



1 <i m <i n 

Proof It is well known the following Taylor series 



(27T)" 



-dm—dn 



(14) 
(15) 



(//2) 
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From and (16), 



Sm,n {f) 



dm H~ 

6 



d rn -d r , 



dn ( f_ 

2 



-dm-d n 



2-{d m +d n ) 



(16) 



2 



Substituting 5 m , n (/) in (12) with (17), we finally get (M4j. 



(17) 



□ 



The asymptotic property of wavelet correlation [T] was also proved for more 
general cases of short memory in |23 1. From (13) and (|14|), the wavelet correla- 



(18) 



tion of x m (t) and x n (t) asymptotically converges to 

Pm,n(j) ~> Pm,n : = D m,nY(d m , d n ) as j Y OO 
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Figure 2: The ratio of fractal connectivity to nonfractal connectivity in a bivari- 
ate ARFIMA(0, <i, 0) process over memory parameters. Note that Hi = d\ + 0.5 
and H 2 = d 2 J r 0.5. 



where 

T(d m ,d n ):=^p^=. (19) 

The asymptotic wavelet correlation p^n * s called the fractal connectivity of 
x m (t) and x n (t). The ratio of fractal connectivity to nonfractal connectivity is 
given from ( p~8| ) by 



D 



T(d m ,d n ). (20) 



The ratio depends just on a pair of long memory parameters as depicted in Fig. 
|2j As the difference of long memory parameters increases, fractal connectivity 
gets away from nonfractal connectivity. On the other hand, fractal connectivity 
is nearly identical to nonfractal connectivity if the difference of two memory 
parameters approaches zero. 

4 Estimation of fractal connectivity and nonfrac- 
tal connectivity 

The multivariate FIP model indicates that both nonfractal connectivity and 
fractal connectivity can be well estimated if memory parameters are known. 
Here we deal with a simple case such that a given time series is approximated 
as a multivariate FIN process (mFIN). The estimation of nonfractal connectivity 
and fractal connectivity is organized as two steps. The first step is to estimate 
memory parameters, and the next step is to compute both fractal connectiv- 
ity and nonfractal connectivity given the memory parameters. In this section, 



several wavelet-based techniques which can be exploited in each step are in- 
troduced. All of these methods are based on the wavelet transform which are 
optimal to investigate the properties of long memory processes. 



4.1 Estimation of memory parameters 
4.1.1 Wavelet least- mean- squares method (LMS) 



By taking logarithm to (14), we obtain the following linear relationship of 
wavelet variance in the logarithm scale: 

log 2 [u m (j)] « 2d m j + Cm. (21) 

It enables us to estimate the memory parameters d m by linear regression over 
a given scale interval J = j\ ow — > jhigh on the basis of the biased estimator of 
wavelet variance given by 

^ = ^1^^) ( 22 ) 

3 t=l 

where rij is the number of coefficients in scale j except boundary coefficients |25 1. 
In a similar manner with |23|, the optimal scale interval J'opt for linear regression 
in (21) is globally determined by 

Jopt = argrmna£ s (J) (23) 
where JJ is the space of all scale intervals, A*/ — Jhigh — Jiow 

+ 1, and 

°ls(J) = ^jEE \ l0 g2 Pm] - 2d m j ~ Cm) • (24) 
j=jiow m=l 

4.1.2 Wavelet maximum-likelihood method (ML) 

The likelihood function for memory parameter d m and asymptotic variance G m 
is given by 

L (d m ,^ m \ Xm (t)) := (27r)JV/2 1 |Sro|1/2 e- T ^ x/2 . (25) 

The matrix S m denotes the covariance matrix of x m (i), and can be replaced 
by S m := W T A>V where W is a wavelet transform matrix and A is a diagonal 
matrix which has diagonal elements A m (j) := Vm(j) given in ( [l4| ) for j = 
1,..., J as an average value of spectral density function (SDF) over the band 
[l/2^' +1 ,l/2^]. The reduced log likelihood function can be obtained on the 
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basis of the Brockwell and Davis's method |26 : 

I (4,7m I Em (*) ) (26) 

= -21ogL^ m ,7 m |x m (t)) -iVlog(27r)-iV 

J 

= TV log (7 m?m ) + log(A m (J + 1)) + 5^- log (A m (j)) 

i=i 

with TV = 2 J , iVy = AT/2^ and 

fm,n 



N VA m ,„(J + l) 
J 1 ^ \ 



where V^, are the scaling coefficients at scale J and W^t is the t-th element 
of j-th level wavelet coefficients. The optimal memory parameter d m can be 



estimated by minimizing (26) with respect to d m |25 1. 



4.2 Estimation of short-memory covariance 
4.2.1 The SDF-based method (SDF) 

The estimator 7 m , n of short memory covariance can be semiparametrically com- 



puted by (27). Since V m ~ in stationary long memory processes, the equation 



can be approximated as 

1 J 1 Nj ~' 
^ W AtEa— m ^ W m (j,t)W n ti,t). (28) 

j=l m ' n ^ ) t=0 

4.2.2 The covariance- based method (COV) 

Alternatively, the short memory covariance can be estimated by exploiting the 
properties such that the sum of wavelet covariances over all scales is identical 
to the covariance of a time series; i.e., 

_ 2 cov(V m ,V n ) cov(W m (j,t), W n (j,i)) 



^2 cuv i^mu,^ yy-nKj^)) ^ 



N ^ 23 

3 = 1 



Since cov (V m , V n ) « for a FIN process, the estimator of short memory covari- 



ance can be obtained from (14) and (29) as follows 



-2 



7m n = ( 2 7T) dM . (30) 
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4.2.3 The linearity-based method (LIN) 

The estimator of short memory covariance can be also obtained in the other 
way based on the linearity of wavelet covariance over scales as follows. 

Jm,n = ~ 2 7„"" 1 , ,j 27T) d - +rf " (31) 

B TO) „cos (I (d m -d n )) 

where 

1 J 

C m , n = J ^2 [log 2 Vm,n{j) ~ (d m + d n )j] , (32) 
J 3 = 1 

-E>m,n := ~Z ~j ~j • (^3) 
1 Cl rn Cl n 

4.3 Estimation of fractal and nonfractal connectivity 

After the estimators for memory parameters d and the short memory covariance 
T are obtained, the nonfractal connectivity D m n can be estimated by using ftlOh 
as follows 

D m , n = ^ . (34) 



Likewise, fractal connectivity n can be estimated from (18) and (34) as fol- 
lows 



5 Simulation study 

In this section, the performance of three wavelet-based estimators for nonfractal 
connectivity was evaluated. We also analyzed the influence of short memory 
condition, dimension, and length of time series on the estimation of nonfractal 
connectivity By combining a connectivity estimator with a memory parameter 
estimator, six pairs of estimator, such as LMS-LIN, LMS-COV, LMS-SDF, ML- 
LIN, ML-COV, ML-SDF methods, were finally tested. 



5.1 Setup 

We simulated multivariate ARFIMA(p, d, 0) processes that belong to the FIP 
model. First, the short memory U(t) in ||5) was given as an ARMA(p, 0) process 
as follows 

U(t) = *; 1 (L)As(t). (36) 



In (36), £i(t) for i = 1, • • • , q is an i.i.d. random variable where 

/nn I 1 if m = n /r . 

cov(e m (t),e n (t)) = \ (37) 
10 llmf n, 
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and 



o 

(i o ••• 

Ola 



(38) 



a 



a 

V 



(39) 



If we set 



1 ± VI - bp 



(40) 



where 6 = p(q— 2) — (g— 3), the short memory correlation is forced to be £> m , n = p 
for m, n > 1 and m ^ n. Afterwards, the memory parameters d were equally 
distributed over d G (—1/2, 1/2), and the multivariate ARMA(p, 0) process was 
filtered by the LM filter defined in 

5.2 Effects of short memory condition 

To study the effects of short memory conditions on the performance of esti- 
mators, we performed Monte Carlo simulations with 100 repetitions of four- 
dimensional ARFIMA(p, 0) processes under four different types of short mem- 
ory condition in (|38|) and (|39|): 



(IA) A = I and ifk,i = 

(IB) A = I, ipk^i — 0.9 and cp kj i = for i > 1 



(2A) A = A and (p k ^ 







(2B) A = A , ifk,i = 0-9 and (p k ,i = for i > 1 



where d = {0.2,0.4,0.6,0.8} and A was set with p = 0.3 in (40). In the 



conditions (1A) and (IB), each short memory process is statistically independent 
of each other while the conditions (2A) and (2B) let short memory processes 
be cross-correlated. On the other hand, the conditions (IB) and (2B) let each 
process be autocorrelated. 

In Fig. [3j all methods were weak in the conditions (IB) and (2B) where short 
memory parts were more auto-correlated; the relative decrease in consistency for 
two cases was common to all six methods. The deteriorated performance in (IB) 
and (2B) is a foreseeable result since these short memory conditions no longer 
follow the assumption of mFIN in (|8| adopted for our proposed estimators. 
Hence, our proposed estimators were not efficient when the set of short memory 
signals cannot be approximated as a multivariate i.i.d. process. 
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(b) ML-LIN 



B 



B 



B * B 



1A 1B 2A 2B 

(c) LMS-COV 



1A 1B 2A 2B 
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(f) ML-SDF 



Figure 3: Box plots of bias in estimation of nonfractal connectivity for simulated 
ARFIMA(p, <i, 0) processes with different short memory conditions. 



In the cases (2A) and (2B) where short memory processes are cross-correlated, 
LMS-LIN, ML-LIN, and ML-COV estimators were not significantly biased. 
However, Fig [4] and [5] show that all estimators except LMS-LIN are more biased 
as the correlation of short memory increases when the dimension (the number 
of time series) is large. 

In these experiments, the short memory processes were given as a multivari- 
ate ARMA(0,0) process with innovation which fulfills (39) and (40). The short 
memory correlation coefficient in (40) was set by either p = 0.2 or p = 0.8. 
In Fig [4] and |5j the absolute bias of estimators except LMS-LIN increased as 
the dimension increases, and the increasing rate of bias was faster when the 
short memory correlation was high. Hence, the high correlation of short mem- 
ory results in the deterioration of estimation performance. All the above results 
show that the bias of estimators tends to be associated with cross-correlation 
of short memory parts, but also that the consistency tends to be related to 
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Figure 4: Box plots of bias in estimation of nonfractal connectivity for simulated 
ARFIMA(0, <i, 0) processes according to variable dimensions when the short 
memory correlation is p = 0.2. 



auto-correlation of each short memory part. In summary, the performance of 
our proposed estimators for nonfractal connectivity is manifestly influenced by 
the short memory conditions. 

5.3 Effects of dimension and length of time series 

Fig [4] and [5] show that the estimators of nonfractal connectivity except LMS-LIN 
depend on the dimension of time series. The LMS-LIN method was relatively 
less biased even in high dimension than other methods, however the consistency 
was large. On the other hand, the other methods were more biased as the di- 
mension increases; especially the increase in bias was more prominent in the 
LMS-SDF and ML-SDF methods. Hence, the increase in the number of time 
series leads to the poor performance of estimating nonfractal connectivity. Nev- 
ertheless, the ML-COV method has the best consistency even in high dimension 
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Figure 5: Box plots of bias in estimation of nonfractal connectivity for simulated 
ARFIMA(0, <i, 0) processes according to variable dimensions when the short 
memory correlation is p = 0.8. 



and high correlation. 

Fig. [6] shows that the performance of nonfractal connectivity estimation is 
associated with the length of time series. In multivariate ARFIMA(0, 0) pro- 
cesses with zero cross-correlation between short memory parts, the consistency 
of all methods was improved as the length of time series increases. However, 
the LMS-LIN method had greater bias than ML-LIN when the length of time 
series was small. In Fig. [3j the LMS-LIN and ML-LIN methods had different 
performance in all tested short memory conditions even though they have the 
common approach to estimating nonfractal connectivity. These results imply 
that the estimation of nonfractal connectivity can be affected by the estimator 
of memory parameter. 

Concluding Fig. [3][6j the ML-COV method would be the best choice as an 
estimator of nonfractal connectivity since it exhibits small bias in various cases 
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Figure 6: Box plots of bias in estimation of nonfractal connectivity for simulated 
ARFIMA(0, <i, 0) processes according to variable length of time series when the 
short memory correlation is p = 0. 



of short memory conditions and high consistency even in high dimension and 
high short memory correlation. As shown in Fig [7| the connectivity matrix es- 
timated by the ML-COV method has smaller difference with the original short 
memory correlation matrix than those produced by other methods have. How- 
ever, it is theoretically expected that the covariance-based method including 
ML-COV may be more sensitive to additive noises than other methods. In the 
case that the signal-to-noise ratio is low, the ML-LIN method may have better 
performance than ML-COV. 
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Figure 7: An example of nonfractal connectivity estimation in a simulated 8- 
dimensional ARFIMA(0, d, 0) process. 

6 Resting state functional MRI 

We applied our proposed estimator of nonfractal connectivity to a resting state 
fMRI data of the anesthetized rat brain taken from 4.7T MRI scanner. We man- 
ually separated the 15 ROIs from the anatomical MRI image, and mapped them 
into the 64 x 64 x 8 volume of blood-oxygen-level-dependent (BOLD) signals in 
fMRI by using the FLIRT (FMRIB's Linear Image Registration Tool). These 
ROIs correspond to aCG, CPu-L, CPu-R, MEnt+MEntV-L, MEnt+MEntV-R, 
HIP-L, HIP-R, Sl-L, Sl-R, S2-L, S2-R, LSI+MS, TE-L, TE-R, and TH [27]. 

We computed both the Pearson correlation and nonfractal connectivity by 
using the ML-LIN method. In Fig. |8j the estimated nonfractal connectivity has 
significantly different patterns from the ordinary correlation. As shown in the 
modularized graph representations of Pearson correlation and nonfractal con- 
nectivity obtained after thresholding the number of edges by 20, the nonfractal 
functional network tends to exhibit increased modularity while the correlation- 
based functional network exhibits high randomness. 
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(a) Nonfractal connectivity 
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Figure 8: An example of nonfractal connectivity estimation in resting state 
fMRI data of the anesthetized rat brain. 

7 Conclusion 

In this article, we modeled a resting state neuroimaging signal as a fraction- 
ally integrated process and introduced the nonfractal connectivity as a novel 
concept of resting state functional connectivity. There is no empirical evidence 
which demonstrates that the nonfractal connectivity reflects the correlation of 
spontaneous neuronal population activities. Through empirical analyses and 
computational modeling of resting state neuroimaging data, the association of 
nonfractal connectivity and neuronal population activities needs to be clarified 
in the future. 

We also proposed several wavelet-based methods for estimating nonfractal 
and fractal connectivity. These estimators are optimal under the assumption 
that the given signals can be approximated by an mFIN while neuroimaging 
signals would have various short memory and can be seriously contaminated by 
physiological or system noises. Hence, the estimators need to be improved in 
a variety of short memory conditions and the existence of additive noises. It 
would be also valuable to figure out the biological mechanism of fractal behavior 
which gives rise to the distortion in functional connectivity. All these challenges 
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may give us insight into the relationship of resting state functional connectivity 
and fractal behavior. 
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